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Abstract 



Using 12.7 fb _1 of data collected with the CLEO detector at CESR, we 
observed two-photon production of the cc states XcO and Xc2 in their decay to 
7r + 7r~7r + 7r~. We measured L 77 (^ c ) xB(x c — * 7r + 7r~7r + 7r~) to be 75 ± 13 (stat) 
± 8 (syst) eV for the XcO and 6.4 ± 1.8 (stat) ± 0.8 (syst) eV for the Xc2, im- 
plying L 77 (x c0 ) = 3.76 ±0.65(stat)± 0.41 (syst) ±1.69(br) keV and L 77 (x c2 ) = 
0.53 ± 0.15(stat) ± 0.06(syst) ± 0.22(br) keV. Also, cancelation of dominant 
experimental and theoretical uncertainties permits a precise comparison of 
r 77 (xco)/r 77 (xc2), evaluated to be 7.4 ± 2.4(stat) ± 0.5(syst) ± 0.9(br), with 
Q CD-based predictions. 
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In this Letter, we report a study of two-photon production of the C-even 1 3 P charmonium 
states Xco and Xc2 using the CLEO detector at the Cornell Electron Storage Ring (CESR). 
The charmonium system is an ideal testing ground for quantum chromodynamics (QCD). 
Perturbative QCD (PQCD) provides predictions for the two-photon widths r 77 (Xco) and 
r 77 (xc2)- These predictions involve charmed quark mass factors, non-perturbative factors, 
and wave-function dependence, all of which cancel in the ratio r 77 (x c0 )/r 77 (x c2 ). The 
experimentally measured quantities are the two products T 77 (x c ) x B(xc —> 7r + 7r~7r + 7r~), 
Xc = Xco or Xc2- Systematic uncertainties in the measurements mostly cancel in the ratio 
of these products. Further, the contribution to the uncertainty on r 77 (x c0 )/T 77 (x c2 ) from 
the uncertainty on the ratio B(xco K + -n~-K + -n~) /B(xc2 /T + 7r~7r + 7r~) is greatly reduced 
because both branching fractions have been measured in the same experiment leading to 
cancelation of their systematic uncertainties in the ratio. Thus the ratio r 77 (x c0 )/r 77 (x c2 ) 
affords a more precise comparison of theory and experiment than do the individual r 77 (x c )- 

The two-photon width of a Xc meson can be determined from a measurement of its two- 
photon cross section. The ratio of the two-photon width to the two-gluon width of a Xc meson 
can be calculated in PQCD with reduced uncertainties due to cancelation of charmed quark 
mass factors, non-perturbative factors, and wave function dependence. In next-to-leading 
order (NLO) PQCD one obtains the following relationships ]IJ : 

r 77 (x c0 ) _ 8a 2 (l + 0.18a s /7r) 

r 99 (Xco) 9«2 (l + 9.5a s /vr) ' 1 ' 



r 77 (Xc2) = 8a 2 (1 - 5.3a 8 /7r) 

T gg (Xc2) 9aUl-2.2a s /nY U 

The width of the Xco meson can be assumed to be dominated by its two-gluon component, 
so T gg (xco) ~ r tot (xco) = 14.9i|J MeV [§]. Using a value of the strong coupling constant 
a s = 0.28 0], one obtains the NLO PQCD prediction r 77 (xco) = 5.0 ± 0.8 keV. Due to 
the uncertainty in the charm mass scale, we also calculate the NLO PQCD prediction at 
a s = 0.35 and find r 77 (xco) = 2.9 ± 0.5 keV. A measurement reported in a thesis gave 
r 77 (Xco) — 4.0 ± 2.8 keV Hj. The E835 collaboration reported an upper limit of r 77 (Xco) < 
3.47 keV (95% C.L.) f|. 

The two-gluon component of the Xc2 width can be extracted from its width r to t(Xc2) = 
2.00 ± 0.18 MeV 0] by subtracting the radiative width T(xc2 — ► an d the color-octet 

width in which the Xc2 decays via three gluons to light hadrons. This latter contribution 
has been shown in Ref. || to be equal to the hadronic width of the xa meson. In the case 
of the Xc2 meson, the color octet contribution to the full width is significant; for the Xco 
meson, it is negligible. Hence, one obtains the NLO PQCD prediction r 77 (xc2) = 0.47±0.04 
keV using a value of a s = 0.28. This prediction becomes r 77 (xc2) = 0.25 ± 0.02 keV for 
a s = 0.35. The Particle Data Group (PDG) value of T 77 (x C 2) is 0.47 ± 0.17 keV @, where 
the uncertainty includes a scale factor of 1.9 to account for the poor consistency among 
the various measurements. The recent measurement of T 77 (Xc2) = 0.270 ± 0.049(stat) ± 
0.033(syst) keV [|J by the E835 collaboration was not included in this average. 

The comparison of the experimental results for T 77 (Xco) and r 77 (x C 2) with the theoretical 
predictions is hampered by the lack of precision in the measured Xc hadronic branching 
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fractions, which are needed to extract the two-photon widths. Exploiting the fact that the 
Xco and the Xc2 mesons were both detected in their decay into 7r + 7r~7r + 7r~, we calculate 
the ratio of their two-photon widths where the dominant branching fraction uncertainties 
cancel. In addition, some of the dominant systematic uncertainties in the r 77 (x c o) and 
r 77 (Xc2) measurements cancel in the ratio because the Xco an d the Xc2 mesons were detected 
in the same experiment. The NLO PQCD prediction for the ratio of the two-photon widths 

is i 

r 77 (Xco) = 15 (l + 0.18a a /7r) 

r 77 ( Xc2 ) 4(1- 5.3a» ' 1 ' 

Here too, charmed quark mass factors have canceled, and we have assumed |^ c0 (0)| 2 = 
|\l/^ c2 (0)| 2 in the non-relativistic limit, where \J r ' (0) denotes the derivative of the Xc wave 
function at the origin. Relativistic corrections due to explicit modification of the meson 
decay amplitudes decrease the above ratio whereas corrections arising from the modification 
of the Xc wave-functions increase the above ratio [[|. In NLO PQCD, predictions for this 
ratio range from 7 to 9 when a s varies from 0.28 to 0.35. Experimentally, this ratio is found 
to be 8.7 ± 6.9, using the single available r 77 (x c0 ) measurement || and the PDG value of 
T 77 (Xc2) 0- Thus a precise measurement of the individual two-photon widths of the Xco 
and the Xc2 mesons and, more importantly, their ratio are of interest for a comparison with 
theory. 

At CESR, which operates at and near the T(4S) resonance, C-even charmonium states 
are produced through the fusion of two space-like photons, radiated by the 5.3 GeV e + and 
e~ beams. The data used in this study correspond to an integrated luminosity of 12.7 ftT 1 
and were collected with two configurations (CLEO II [0 and CLEO II. V |j) of the CLEO 
detector. Approximately one third of the data was taken with the CLEO II configuration. 
The detector components most useful for this study were the concentric tracking devices 
for charged particles, operating in a 1.5 T magnetic field generated by a superconducting 
solenoid. For CLEO II, this tracking system consisted of a 6-layer straw tube chamber, 
a 10-layer precision drift chamber, and a 51-layer main drift chamber. The main drift 
chamber also provided measurements of the specific ionization loss, dE/dx, used for particle 
identification. For CLEO II. V, the straw tube chamber was replaced by a 3-layer, double- 
sided silicon vertex detector and the gas in the main drift chamber was changed from a 
50:50 mixture of argon-ethane to a 60:40 helium-propane mixture. These changes gave rise 
to a significant improvement in the momentum and dE/dx resolutions for charged tracks. 
Photons were detected using the high-resolution electromagnetic calorimeter consisting of 
7800 Csl crystals. The Monte Carlo simulation of the CLEO detector response was based 
upon GEANT ||. Simulated events were processed in the same manner as the data to 
determine the Xc —> 7r + 7r~7r + 7r~ detection efficiencies and the 7r + 7r~7r + 7r~ mass resolutions 
at the two Xc meson masses. 

In the two-photon process e + e~ — > e + e~77 — > e + e~Xc e + e~7r + 7r~7r + 7r _ , the photon 
propagators dictate that the photons are almost real ("on shell"). Therefore the incident 
leptons are scattered at very small angles with the beam and remain undetected. Such 
"untagged" events typically have low transverse momentum {px = | J2i Pr* | , i = 1 — 4) and 
low visible energy. 
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The events were required to have exactly four charged tracks with zero total charge. The 
background from processes other than two-photon production was suppressed by requiring 
that the \c candidate reconstructed from these four charged tracks has transverse momentum 
less than 0.4 GeV/c and that the visible energy in the event be less than 6.0 GeV. A x 2 
probability was constructed using the dE/dx information from the four tracks, and required 
to be greater than 10%. Also, because the final state had no expected energy deposits in 
the calorimeter from neutral particles, the total calorimeter energy in the event not matched 
to charged tracks was required to be less than 0.6 GeV. We required that at least two of 
the four tracks traversed all layers of the tracking volume to achieve a well-modeled trigger 
simulation in the Monte Carlo. 




3200 3400 3600 3800 

M (4jr) (MeV) 

FIG. 1. The 7T + 7r - 7r + 7r~ invariant mass (data point with errors). The solid line is the fit with 
a xVd.o.f. = 44/54. 

The invariant mass distribution of the selected four-pion events is shown in Figure |I[ 
Enhancements near the known Xco and Xc2 masses are seen upon a smooth background. 
To obtain the number of signal events and mass values, the mass distribution was fitted 
as follows. The background was fitted with a power law function: A ■ W™, with W 77 the 
7r + 7r~7r + 7r~ invariant mass and A and n parameters to be fitted. The Xco signal was fitted 
with a spin-0 relativistic Breit-Wigner function with a width equal to the PDG value of 
14.9 MeV (describing the natural line shape) convolved with a double Gaussian function 
(describing the detector resolution). The Xc2 signal was fitted using only a double Gaussian 
resolution because the PDG value of 2.0 MeV for its width is negligible in comparison to the 
mass resolution. The parameters of the resolution functions were determined from Monte 
Carlo samples of the Xco and Xc2, generated with zero width. The invariant mass resolution 
was approximately 9 MeV for both the resonances. We performed a simultaneous, binned, 
maximum-likelihood fit to the invariant mass distribution using these three functions. 

From the fit, we obtained 234 ± 40 candidate Xco events and 89 ± 25 candidate Xc2 
events, with masses of 3416.5 ± 3.0 MeV and 3559.9 ± 2.9 MeV, respectively. Our measured 
masses are consistent with the PDG values, 3415.0 ± 0.8 MeV and 3556.18 ± 0.13 MeV @ 
respectively. 

Comparison of the two-photon crosssection of a resonance measured from data with that 
estimated from Monte Carlo for a given two-photon width allows extraction of the r 77 (x c ). 
The Monte Carlo crosssection (<tmc) estimation was based upon the BGMS formalism [|TUJ. 
The cross-section in data was calculated by dividing the number of events from the fit 
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TABLE I. Systematic uncertainties in the two T~y(xc) x &(Xc ~^ vr + 7r _ 7r + 7r _ ) measurements. 
The total systematic uncertainties were obtained by adding the individual contributions in quadra- 
ture. Most systematic uncertainties are pairwise correlated and are marked with an asterisk. 



Source of Uncertainty 


XcO (%) 


Xc2 (%) 


Event pt* 


7 


7 


Particle Identification* 


4 


6 


Unmatched Neutral Energy* 


4 


4 


Helicity 





4 


Mass Calibration* 


1 


3 


Width of Resonance 


4 





p° Substructure* 


2 


2 


Trigger* 


2 


2 


Tracking* 


2 


2 


Detector Resolution* 


1 


1 


Signal Shape Parameter* 


2 


2 


o"mc (inc pole mass)* 


1 


1 


Luminosity* 


1 


1 


Total 


11 
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by the detection efficiency, the integrated luminosity, and the branching fraction B(xc ~ > 
7r + 7r~7r + 7r~). The detection efficiencies were 19.6% and 21.7% for the Xco and the Xc2, 
respectively. Due to the large uncertainties in the branching fractions, we preferred to 
express our primary results as the products of the branching fraction and the two-photon 
width. We obtained r 77 (x c0 ) x B(xco 7r + 7r~7r + 7r~) = 75 ± 13(stat) ± 8(syst) eV and 
T 77 (Xc2) x £>(Xc2 — ► 7r + 7r~7r + 7r~) = 6.4 ± 1.8(stat) ±0.8(syst) eV. From these primary results 
we obtained the ratio r 77 (x c0 )/r 77 (x c2 ) = 7.4±2.4(stat) ±0.5(syst) ±0.9(br), where the last 
uncertainty corresponds to branching fraction uncertainties that do not cancel in the ratio. 
Using the known branching fractions, B(xco ~^ 7r + 7r~7r + 7r~) = 2.0 ± 0.9 % and B(xc2 
7r + 7r~7r + 7r~) = 1.2 ± 0.5 % 0, we obtained the two-photon widths as r 77 (x c0 ) = 3.76 ± 
0.65(stat)±0.41(syst)±1.69(br) keV and r 77 (x c2 ) = 0.53±0.15(stat) ±0.06(syst) ±0.22(br) 
keV, where the last uncertainties correspond to the systematic uncertainties arising from the 
hadronic branching fraction. 

Sources of systematic uncertainty in the measured product branching fractions are sum- 
marized in Table I. The uncertainties which were correlated are marked by an asterisk. The 
uncertainties were added in quadrature to obtain the total systematic uncertainty. The sys- 
tematic uncertainty of the ratio of the two-photon widths takes into account correlations 
among the uncertainties. Systematic uncertainties in this ratio were separately estimated 
for our experimental measurement and the branching fraction measurements from other ex- 
periments pT|-4T3]]. Though the individual measurements of B(x c o — * 7r + 7r - 7r + 7r~) and the 
B{Xc2 — * 7r + 7r~7r + 7r~) have poor consistency, the ratios of these two branching fractions are 
consistent. 
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The systematic uncertainties were dominated by the uncertainties in the simulation of 
the signal pt distribution, the particle identification procedure and the simulation of the 
unmatched neutral energy. 

In our cross-section estimate from Monte Carlo we had assumed that the Xc2 meson was 
produced only in the helicity 2 state. The helicity state has zero two-photon width in 
the non-relativistic approximation, while in relativistic models it is predicted to be up to 
4% of the helicity 2 state ||. We allowed the mass of each resonance to vary between our 
measurement and the PDG value to estimate the systematic uncertainty from uncertainty 
in the mass scale. We similarly varied the width of the Xco within its known uncertainty. 
Using the known branching fractions of the Xco an d Xc2 [0 into p°tt + 7c~ and estimating the 
difference in detection efficiencies using Monte Carlo we assigned the systematic uncertainty 
due to the presence of resonant p° substructure in the Xc decay. Ref. [O had reported no 
significant p°p° sub-structure in the Xco and Xc2 decay. 

We investigated the possible effects of interference between the four-pion decay of the Xc 
states and non-resonant production of four pions in two-photon interaction. From Monte 
Carlo simulation we found that a fifth of the selected events were from annihilation r-pair 
production. Another one-third were estimated to be not of two-photon origin based on the 
candidate transverse momentum distribution. Neither of these two classes of events can 
interfere with signal events. Interference can distort the signal shape derived from the signal 
Monte Carlo. We observed no such distortion demonstrated by the good fit x 2 m t ne two 
regions of the resonances: 12.4 for 12 degrees of freedom. 

The systematic uncertainty of 0.5 in the ratio r 77 (x c0 )/r 77 (x c2 ) takes the correlations 
marked by an asterisk in Table I into account. The systematic uncertainty of 0.9 from the 
ratio of branching fractions was calculated taking into account correlations as ascertained 
from Ref. [|ll|] and [ 12 1 . The ratio of the branching fractions has significantly lower relative 
uncertainty than the branching fractions themselves due to partial cancelation of correlated 
systematic uncertainties. The branching fractions B(ip(2S) — *> 7Xc) needed were taken from 



Ref. ||13|| , taking into account correlations between them. 

Our result r 77 (x c0 )/r 77 (x c2 ) = 7.4±2.4(stat) ±0.5(syst) ±0.9(br) represents a significant 
improvement upon the current value of 8.7 ± 6.9 and is compared with the NLO PQCD 
prediction (Eq. 3) in Figure [|. There is good agreement for reasonable values of a s ; however 
one wishes for increased theoretical and experimental precision. Our r 77 (x c0 ) measurement 
of 3.76 ± 0.65 (stat) ± 0.41 (syst) ± 1.69 (br) keV is a significant improvement in experimen- 
tal precision over the single available measurement of 4.0 ± 2.8 keV ||, but limited by poor 
knowledge of the branching fraction. The NLO PQCD prediction of 3 - 5 keV is consistent 
with our measurement. Our r 77 (x c2 ) measurement of 0.53 ± 0.15 (stat) ± 0.06 (syst) ± 
0.22 (br) keV is consistent with the PDG value of 0.47 ± 0.17 keV [§. The PDG average of 
r 77 (Xc2) neglects correlated systematic uncertainties between experiments. The NLO PQCD 
prediction of 0.25 - 0.47 keV is consistent with our measurement and the PDG value. Our 
measurements show that PQCD based calculations are able to predict the ratios of decay 
rates, where non-perturbative effects cancel; however the uncertainties in both theory and 
experiment limit the precision of such comparisons. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Research Corporation, the Natural Sciences and 
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FIG. 2. Comparison of the measurement of the ratio r 77 (x c o)/r 77 (Xc2) (shaded region corre- 
sponds to ± 1 cr) with its prediction from NLO PQCD as function of a s . 
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